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Abstract 


Tropical Rainfall Measuring Mission (TRMM) Microwave Imager (TMI) 
radiometer brightness temperature data in the 85 GHz channel (T85) reveal distinct 
local minima (T85 m i n ) in a regional map containing a Mesoscale Convective System 
(MCS). This is because of relatively small footprint size (~ 5.5 km) and strong 
extinction properties in this channel of the TMI. A map of surface rain rate for that 
region, deduced from simultaneous measurements made by the Precipitation Radar 
(PR) on board the TRMM satellite, reveals that these T85 m in/ produced by scattering, 
correspond to local PR rain maxima. Utilizing the PR rain rate map as a guide, we 
infer empirically from TMI data the presence of three different kinds of 
thunderstorms or Cbs. These Cbs are classified as young, mature, and decaying 
types, and are assumed to have a scale of about 20 km on the average. Two 
parameters are used to classify these three kinds of Cbs based on the T85 data: a) the 
magnitude of scattering depression deduced from local T85 m in an d b) the mean 
horizontal gradient of T85 around such minima. Knowing the category of a given 
Cb, we can estimate the rain rate associated with it. Such estimation is done with 
the help of relationships linking T85 m i n to rain rate in each Cb type. Similarly, a 
weak background rain rate in all the areas where T85 is less than 260 K is deduced 
with another relationship linking T85 to rain rate. In our rain retrieval model, this 
background rain constitutes stratiform rain where the Cbs are absent. Initially, these 
relationships are optimized or tuned utilizing the PR and TMI data of a few MCS 
events. After such tuning, the model is applied to independent MCS cases. The 
areal distribution of light (1-10 mmhr'l), moderate (10-20 mmhr'l), and intense 
(> 20 mmhr 1 ) rain rates are retrieved satisfactorily. Accuracy in the estimates of the 
light, moderate, and intense rain areas and the mean rain rates associated with such 
areas in these independent MCS cases is on the average about 15 %. Taking 



advantage of this ability of our retrieval method, one could derive the latent heat 
input into the atmosphere over the 760 km wide swath of the TMI radiometer in 
the tropics. 


- *, 


1. Introduction 

Several microwave radiative-transfer rain retrieval models (see for e.g. Wu 
and Weinman, 1984; Kummerow et al., 1989; Smith and Mugnai, 1992; and 
Prabhakara et al., 1995) have taken into consideration microwave scattering by dense 
solid ice particles, which are generally associated with convective clouds. However, 
Schols et al. (1999) indicate that in order to explain satellite-borne Special Sensor 
Microwave Imager (SSM/I) 1 brightness temperature data associated with stratiform 
precipitation produced by nimbostratus clouds, radiative properties of frozen and 
melting ice aggregates (snow) with low density have to be considered. Otherwise, a 
significant uncertainty in estimation of such stratiform rain rates from SSM/I data 
will result. In order to objectively sense dense ice particles (frozen rain drops and 
graupel) or less dense ice particles (frozen and melting ice aggregates), a spectral 
signature is needed from the microwave radiometer data. Such a distinct signature 
is missing in the limited information contained in the dual-polarization, 
multispectral SSM/I data (Schols et al., 1995; Prabhakara et al., 1998). This is because 
measurements in all the channels of this instrument have a high degree of 
redundancy. This lack of information in the microwave radiometer data can lead to 
significant errors in rain retrievals that are based solely on radiative transfer theory. 


1 Radiometer of the Defense Meteorological Satellite Program (DMSP) that has 19, 37 and 85 GHz 
channels in dual polarization and a 22 GHz channel in vertical polarization. This radiometer observes 
the earth's surface and atmosphere in a conical scan with an incidence angle of ~ 50° (for more details 
see Hollinger et al, 1985 ). 
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This problem associated with retrieving rain using microwave radiative 
transfer models is exemplified over oceanic regions by the studies of Heymsfield et 
al. (1996) and McGaughey and Zipser (1996). They find from aircraft radar and 
microwave radiometer observations that retrievals from current models that 
consider dense ice particles perform poorly. Based on the Schols et al. (1999) study, 
we infer that this is largely due to improper specification of the nature of the ice 
particles in the models. We may remark that the spatial resolution of these aircraft 
instruments is a few kilometers, which is much smaller compared to that of the 
SSM/I. Crude spatial resolution of the satellite-borne microwave radiometers is 
often cited as the dominant reason for significant errors in rain retrievals (see for 
e.g. Kummerow, 1998). However, we contend, based on the errors in rain rate 
retrievals from SSM/I and aircraft radiometer data, that spatial resolution is not the 
only cause of such errors. The ability to discriminate the nature of ice hydrometeors 
is also an important factor. 

On land regions, rain retrieval from microwave radiometer data is further 
complicated by variable surface emissivity introduced by vegetation, surface 
wetness, water bodies and terrain. Conner and Petty (1998) indicate that the spectral 
information contained in the SSM/I data over land can explain only about 30% of 
the variance present in the radar or the rain gauge observations over land. From 
the above discussion, we conclude that because of the spatial inhomogeneities and 
the uncertainties in the nature of the hydrometeors, radiative transfer models do 
not succeed in retrieving rain rate over land. Thus, in order to get good rain 
retrievals over land and ocean, more information in addition to the microwave 
radiometer spectral measurements is needed. 
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In the studies of Prabhakara et al. (1998, 1999), the rain area in a mesoscale 
geographic region (~ 2° lat x 3° Ion) deduced from the SSM/I data is provided as 
additional information to estimate the mesoscale area-average rain rate. Those 
studies indicate that the rain area in a mesoscale region serves as a more important 
factor in estimating area-average rain rate than the spectral information given by 
the SSM/I data. This mesoscale area-average rain retrieval approach is widely used 
and has its basis in the Area Time Integral concepts of Doneaud et al. (1984) and 
Lopez et al. (1989). A useful inference from these studies is that the atmospheric 
dynamics that lead to the total volume of rain in a mesoscale region are coupled to 
the rain area in that region. 

Based on a similar conceptual framework, in this study the spatial, and the 
spectral, information given by the Tropical Rainfall Measuring Mission (TRMM 2 ) 
Microwave Imager (TM1) radiometer data are used in an empirical fashion to 
identify the presence of thunderstorms (Cbs) and obtain convective rain rate in the 
limited area of those Cbs. The rain areas outside of such limited convective regions 
are considered to be stratiform in nature. The empirical method developed in this 
study is designed to perform well over both land and ocean. Such spatial 
information, derived from satellite infrared radiometer observations, has been 
utilized empirically in an earlier study by Adler and Negri (1988) for the purpose of 
obtaining convective and stratiform rain. However, since microwave data represent 
properties of the dense layers in precipitating clouds, and not just the thin layer near 


2 Tropical Rainfall Measuring Mission. On board the TRMM satellite, there is the cross-track-scanning 
Precipitation Radar (PR) and the conical-scanning TRMM, Microwave Imager (TMI) (see Simpson et al, 
1996 ). 


3 


the cloud tops, microwave techniques perform better than infrared techniques in 
retrieving instantaneous estimates of rain rate (Ebert et al., 1996). 

Hong et al. (1999) developed an empirical rain retrieval method that utilizes 
the spatial distribution of the SSM/I microwave radiometer measurements on 
ocean. This method was optimized with the help of the TOGA-COARE radar 
observations. In this formulation, convective/stratiform index (CSI) for each pixel 
is given by a combination of two quantities: 1) variability in brightness temperature 
near the pixel, and 2) the contrast in brightness temperature at the pixel relative to 
clear sky conditions. This CSI scheme of Hong et al. (1999) has difficulty to limit 
convective rain area to the thunderstorm scale. This scheme is used in the TRMM 
Project as a standard procedure to retrieve rain rates on ocean from the TMI 
microwave radiometer measurements. Olson et al. (1999) have considered in detail 
the merits of the Hong et al. (1999) rain retrievals. They find that this algorithm 
leads to a high bias in rain rate estimation with respect to radar, because it over 
estimates rain in the low to moderate range (0.1 to 16 mmhr 1 ). Hong et al. (1999) do 
not demonstrate the utility of their rain retrieval formulation on land regions. The 
objective of this study is to improve the state-of-the-art, such that we can get 
quantitative convective and stratiform rain retrievals over both ocean and land. 

In the TRMM mission, the TMI measurements cover a swath width of 
760 km. On the other hand, the Precipitation Radar (PR) on board the TRMM 
satellite has a cross-track scan that covers only the 220 km central portion of the TMI 
swath (Simpson et al., 1996). This study aims to develop a method than can exploit 
this much larger scanning capability, and thereby expand temporal and spatial 
coverage of rain information over land and ocean regions beyond that of the PR. 


2. Model to Retrieve Rain Rate Maps with Convective and Stratiform 
Discrimination Based on TMI data 

(a) General Considerations 

The TMI radiometer is similar to the SSM/I, except that it has two additional 
spectral channels (vertical and horizontal polarization) at 10 GHz. Furthermore, 
because of the lower altitude of the TRMM satellite, the spatial resolution of the 
TMI radiometer is about two times finer compared to SSM/I. This finer resolution, 
in conjunction with the spectral information, offers an opportunity to discriminate 
convective rain from stratiform rain. 

We may remark that the vertical resolving capability of the radar gives it a 
distinct advantage over the passive microwave radiometer. Furthermore, radar 
back scatter measurements have a d 6 dependence on the rain drop diameter (Battan, 
1973), while radiometer observations tend to have about a d 3 dependence (see for 
e.g. Olson, 1996). Both of these properties of the radar help it to isolate the rain layer 
and to get a more robust estimate of rain rate. However, it may be remarked that 
rain rate estimates based on radar data are linked to the rain drop size distribution 
(DSD). 


There is a difference in DSD between convective and stratiform rain. Tokay 
and Short (1996) show with disdrometer data that the mode radius of the DSD for a 
given rain rate in stratiform rain is larger. This difference, if it is not accounted for, 
could lead to an overestimation of stratiform rain rate by about a factor of two. 
However, difference in DSD associated with convective and stratiform rain can be 
inferred with the help of horizontal inhomogeneity in rain rate (Steiner et al., 1995). 
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From such inference, errors in rain rate estimation can be minimized. In addition 
to this problem, there could be other temporal variabilities in DSD. Probable error 
in rain rate due to such DSD variability is estimated to be within ~ 50 %. This error 
estimate is based on simulations of radar backscatter using disdrometer observed 
DSD (R. Meneghini, personal communication). Because of these non-negligible 
errors, radar rain rates are calibrated with respect to surface rain gauges. An example 
of such calibration is that performed over Japan with the Automatic Meteorological 
Data Acquisition System (see for e.g. Oki et al., 1997). 

Unlike radar observations, spectral information from measurements made by 
microwave radiometers have poor vertical resolution and ability to discriminate 
stratiform versus convective rain. Therefore, as shown in the study of Heymsfield 
et al. (1996), there are much larger errors associated with rain rates retrieved using 
radiometer spectral data. For these reasons, until better ground truth is available, 
we will use the TRMM Precipitation Radar (PR) rain estimates as a basis to develop 
our TMI rain retrieval method. When better ground truth is available, this TMI 
algorithm can be readily tuned with the help of such information. 

We may point out that the current rain retrieval algorithm used to obtain the 
TRMM Project TMI Standard Rain Rate, RtmI/ on ocean is based on an empirical 
scheme of Hong et al. (1999). This empirical scheme is developed by a tuning 
procedure using SSM/I and radar observations made over TOGA-COARE. This 
procedure contains information from the 19 and 37 GHz channels of the radiometer, 
which are contaminated by surface emissivity effects on land. We contend, for this 
reason, that the scheme of Hong et al. is not adaptable to land regions. With the 
help of the radar data, in this study, we wish to develop one empirical rain retrieval 
scheme that is applicable to both land and ocean. 
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(b) Model Development 


Studies of Houze and Betts (1981) and Houze (1993, 1997) show that in 
Mesoscale Convective Systems (MCSs) there are thunderstorms or cumulonimbus 
clouds (Cbs) in different stages of growth. A Cb may contain a single cloud cell or 
multiple cloud cells. The scale of these Cbs may vary somewhat, but on average it is 
about 20 km. In a gross fashion, a given Cb may be categorized as young, mature, or 
decaying according to its age. The young Cbs have cloud cells with moderate to 
strong vertical motions of several ms -1 , in which rain drops grow primarily through 
the coalescence mechanism. Because these young Cbs have not grown deep into the 
atmosphere, there are relatively small amounts of ice particles in these Cbs. The 
mature Cbs have strong vertical motions, which reach high altitudes and produce 
larger amounts of frozen rain drops and graupel. In the studies of Houze, moderate 
and intense rain rates are reflected by the strength of the narrow updrafts in the 
young and mature Cbs, respectively. For both the young and mature Cbs, the 
strength of the vertical motions and the rain rate are found to decrease rapidly 
toward the Cb edge. On the contrary, it is found in the decaying Cbs that the vertical 
motions are relatively weak (< 1 ms -1 ), which allows sufficient time for the growth 
of falling ice crystals by deposition of water vapor. In turn, these ice crystals form ice 
aggregates before entering the melting layer, and then melt to form rain below this 
layer. This type of rain is stratiform in nature and relatively weak, but is more 
homogeneous in a given area. We take advantage of these concepts of cloud 
dynamics and microphysics in our empirical model, even though we do not have 
observations of vertical velocities and hydrometeors. The model adopted in this 
study to classify Cbs into the three categories, as described below, depends on T85 and 
its spatial distribution. 


In Figure la, a map of a rain event over the Southeast United States on 
5 June 1998 (Day 156) is depicted using the TMI T85 horizontal polarization data. 
The horizontal polarization is preferred because of the surface emissivity effect, 
which gives it a larger dynamic range over land and ocean compared to the vertical 
polarization data. We emphasize the 85 GHz TMI data in this study because, of all 
the channels of this radiometer, the 85 GHz channel has the finest spatial resolution 
and the highest extinction to ice and liquid hydrometeors. In particular, the 
scattering of microwave radiation, which according to radiative transfer theory 
results in a depression in the 85 GHz data, is strongly related to cloud ice content 
(Wu and Weinman, 1984). The surface rain rate map deduced from the PR data for 
the same time period (within two minutes) is shown in Figure lb. These maps 
represent the 220 km central portion of the TMI swath. 

The PR surface rain rate map shown in Fig. lb has a character similar to rain 
maps derived from conventional radars (see for example Houze, 1993). In 
particular, the striking contrast between heavy rain rate in convective cells and 
weak rain rate present in the extensive stratiform region is faithfully reproduced in 
the PR data. The commonly observed log-normal type of frequency distribution of 
rain rate (Jones and Sims, 1978) is also revealed by the PR surface rain rates. This 
frequency distribution is illustrated in Figure lc using the PR surface rain rates 
displayed in Fig. lb. These properties of the PR data reinforce their credibility. 

From the map of PR surface rain rate shown in Fig. lb, we discern that there 
is a squall line composed of several intense Cbs. Furthermore, there are several 
other weaker Cbs. Some of these weaker Cbs are present in the extensive stratiform 
region. By comparing the PR rain rate map with the map of T85, we notice there are 
local minima in T85 (T85 m in), which correspond roughly to rain maxima in the PR 
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rain map. These rain maxima represent the thunderstorms or Cbs. We may point 
out that a local minimum in T85 can be defined in the TMI data with respect to the 
brightness temperature values of its four neighboring footprints given by the TMI 
conical scans, as shown in Figure 2. Apparently, T85 m in can be produced in the 
radiometer footprints by dense ice particles that can strongly scatter microwave 
radiation in the narrow convective cells; or they can be produced by frozen and 
melting ice aggregates that can weakly scatter microwave radiation over a broader 
area. We note from Fig. la that we cannot identify the three different types of Cbs 
based simply on the magnitude of T85 m j n . In order to do such identification, we 
need to infer the spatial structure of the hydrometeors in the Cb. This structure can 
be assessed from the mean of the absolute values of the spatial gradients, i.e., 
| dT85/dr |, within each Cb. Here, r represents the distance between a point in the Cb 
and the location of T85 m j n . This | d T85/ dr | is determined utilizing the four 
neighboring footprint values as illustrated in Fig. 2, and it is analogous to the 
Laplacian (V T85) with an additional condition that T85 at the center point is a 
minimum. 


Given the information T85 m in and |dT85/dr|, we have developed, using PR 
rain rate as a guide, the following simple empirical model to identify young, 
mature, and decaying Cbs. First, because of narrow and relatively strong updrafts in 
the cloud cells of young and mature Cbs that contain dense ice hydrometeors, it is 
assumed that the mean gradient around T85 m i n , |dT85/dr|, is large compared to that 
in the decaying Cbs. Thus, in this model we use a criterion based on | dT85/dr| to 
separate decaying Cbs from convectively active young and mature Cbs. Second, as a 
result of larger amounts of frozen rain drops and graupel, we expect mature Cbs to 
have magnitudes of T85 m i n that are small in comparison to those found in young 
Cbs. Hence, the magnitude of a T85 minimum is used to separate the young and 
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When TQ<Ti,T 2 ,T3and T 4 
Then T m j n =TQ 



A(T*| +T 2 ) + BO3 + T4) - 4 Tq 

’ ” 4 


Figure 2: Illustration of the method to determine the location of a T85 minimum, T85 m in, and the mean spatial 

gradient, | dT85/dr | , around this minimum. The solid curves represent the conical scan of the TMI radiometer, 
and the circles represent the footprints. Along the scan, the separation between footprints is about 4.6 km, while 
the separation between scans is about 13.9 km. In the illustration, the coefficients A and B correspond to 
appropriate weighting factors. 



mature Cbs. The empirical scheme adopted here to separate the three types of Cbs is 
illustrated in Figure 3. 

In Fig. 3, a scatterplot of |dT85/dr| versus T85 m i n is shown. It is based on 
several Cbs from the sample of twelve dependent mesoscale rain events listed in 
Tables 2a and 3a. Six of these events are on land, and the other six are on ocean. As 
shown in this figure, we are assuming that the line ab representing |dT85/dr| 

equal to 1 Kkm- 1 separates the decaying Cbs from the young and mature 
convectively active Cbs. The decaying Cbs are assumed to have |dT85/dr| less than 
1 Kkm -1 , while the active Cbs have | dT85/dr | greater than or equal to 1 Kkm' 1 . The 
simple reason for this assumption is that we expect the ice hydrometeors will be 
relatively homogeneous in space in the decaying Cbs. On the other hand, the 
convectively active Cbs are not expected to have such homogeneity. From the 
figure, we may note that 1 35 m i n is generally greater than 190 K for the decaying Cbs. 
This is not true of the convectively active Cbs. The young and mature Cbs are 
partitioned utilizing the threshold T85 m i n equal to 210 K. This threshold is shown 
in Fig. 3 with the vertical line ef . However, we may point out that the 210 K value 
of this threshold is not critical. We find that within a reasonable range of the 
threshold (210 K > T85 m m > 190 K), our algorithm can be adjusted to get satisfactory 
rain retrievals. These criteria adopted to discriminate young, mature, and decaying 
Cbs are in broad agreement with radiative transfer models that consider dense ice 
particles and less dense ice aggregates. 

In order to illustrate the correspondence between the Cbs indicated by the PR 
and the T85 data respectively shown in Figs, la and lb, a map of the locations of 
T85 min for the rain event of 5 June 1998 (Day 156) is shown in Figure Id. The 
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magnitudes of T85 m i n and | dT85/dr |, and the category of each Cb, are also shown in 
Fig. Id. 

We may note that, because factors such as evaporation below cloud, cloud 
liquid, and fractional cloud cover are highly variable, one cannot determine 
objectively a radiometer measured value of T85 which corresponds to a condition of 
zero rain at the surface. In our rain retrieval model, we assume that surface rain is 
zero in areas where T85 is greater than 260 K. Furthermore, all of the Cbs are 
assumed to have values of T85 m in less than 255 K. This is done to avoid weak Cbs 
that have T85 m i n between 255 K and 260 K that lead to noise in the rain retrievals. 
Note also that in the present study the polarization in the 85 GHz, i.e., P85, is used to 
discriminate open ocean and wet surfaces on land from precipitating clouds where 
T85 is less than 260 K. A value of P85 greater than 15 K is used to screen out these 
surfaces. All rain rates discussed in this study (Rtmi - TMI Standard Rain Rate; Rpr - 
PR Standard Rain Rate; and R - Rain Rate derived from the present scheme) refer to 
rain rate near the surface. 

In the present rain retrieval model, in order to blend the stratiform and 
convective rain rates in a smooth fashion, a weak background rain rate, Rs, is 
assumed to exist in all the areas where T85 is less than 260 K. This background rain 
is considered to be of stratiform type in the areas where Cbs are absent. For this 
purpose, the following linear relationship between the background rain rate, Rs, 
and T85 is used: 

R S = (260 - T85) {dR PR /dT85}. (1) 

Based on radar and radiometer observations, the sensitivity factor dR PR /dT85 is set 
equal to 0.12. 



As stated above, the rain rate in each Cb is superimposed on the background 
rain rate Rg. Each Cb is assumed to have a circular area with a radius of 10 km, with 
its center coinciding with the location of the pixel that has the minimum value 
T85 m in- Note, one T85 pixel that defines a local minimum has an area of about 
30 km 2 . This local minimum leads us to infer the presence of a Cb extending over 
an area of about 310 km 2 , which is about 10 pixels in our model. So the influence of 
a convective Cb is not limited to one pixel. Any pixel within the Cb is assigned an 
appropriate convective rain rate. Pixels outside of a given Cb are considered to 
have stratiform rain regardless of their T85 value. With the 5-pixel scheme utilized 
to identify local minima (see Fig. 2), we can find along a scan line two Cbs with one 
common pixel between them. Thus, there can be two Cbs along a radiometer scan 
line with a center-to-center separation distance of 9.2 km. In such a case, the two 
Cbs, each of which has a 10 km radius, will be overlapping. When there are such 
Cbs with different rain intensities overlapping with one another, we assume the 
rain rate in the overlapping area is governed by the more intense Cb. Thus, there 
can be several overlapping Cbs in a rain scene that can form large areas of 
convective rain. 


We determine the Cb type as described earlier (see Fig. 3). Then based on the 
type, we estimate the mean rain rate in a Cb, excluding the background rain rate Rg, 
using one of the following three linear relationships: 


Young Cbs 

R Y = Fio [(255 - TBS™, ) { dR PR /dT85 min h l 

(2a) 

Mature Cbs 

R m = FlO [ Ry + (210 - T85 min ){ dRp R /d T85 min } M ] 

(2b) 


where Ry = (255 - 210 ) { dRp R /dT85 min }y; and 


Decaying Cbs 

RD = Fio [(255 - T85min) ( dRp R /dT85 m j n }dJ* 

(2c) 
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where RY/ R^/ and Rp) stand for mean rain rates of young, mature and decaying 
Cbs, respectively. This rain rate is superimposed on R$ to get total rain rate in each 
Cb. In Eqs. 2a and 2b, the sensitivity factors (dRpR/dT85 m j n )Y and 
(dRpR/clT85 m j n )M correspond to the convectively active young and mature Cbs, 
respectively. The term Rpr in these equations represents the radar derived mean 
rain rate over the circular area of the Cb, which has a radius of 10 km. As shown in 
Eq. 2c, another sensitivity factor (dRpR/dT85 m j n )D is used for the decaying Cbs. 
These sensitivity factors are optimized with the help of the TMI and PR 
observations of a small sample of MCS cases. Over the ocean, the term Fio is set to 
vary linearly between 0.0 and 1.0 as the 10 GHz brightness temperature varies 
between 100 and 200 K, respectively. This can be expressed as 


F10 = < 


0.0 

(T10-100)/100 

10 


T10 < 100K 
100K < T10 < 200K . 
T10 > 200K 


(2d) 


The term Fio in Equation 2d is assigned a value of 1.0 over land regions, since T10 is 
greater than 200 K. In Table 1, we are presenting the magnitude of dRpR/dT85 min 
for all Cb types. 


As already indicated, the present rain retrieval model based on TMI data can 
estimate the mean rain rate in the Cbs on a scale of 20 km. In order to obtain the 
pattern of rain rate R everywhere in an MCS, including inside the Cbs, we use a 
simple interpolation procedure. As an example, we show below the interpolation 
equation applicable to young Cbs: 

R =R s + [Ry-{ RY/(T85 Ma x-T85)} (T85-T85)]. (3) 

Note that this interpolation scheme does not alter the average rain rate Ry over the 
Cb area. Similar equations are applied to the mature and decaying Cbs. Here, T85 


13 


stands for the 85 GHz brightness temperature at any given point inside a Cb, and 
T85Max and T85 stand for the maximum and mean values of T85 in that Cb. 
Outside a Cb, we set R = R$. Utilizing Eqs. 1, 2 and 3, in the present method we 
have derived the rain rate R everywhere in an MCS where T85 is less than 260 K. 

We find the relationship of Rpr versus T85 m j n for convective young and 
mature Cbs is significantly nonlinear over a wide 85 GHz brightness temperature 
range of about 255 K to 160 K. For simplicity, we fit this nonlinear relationship with 
two linear segments over this range. These two segments are separated at T85 equal 
to 210 K. We find for the decaying Cbs that the Rpr versus T85 m j n relationship is 
weakly nonlinear, and so we represent it with one linear approximation over the 
entire temperature range. The sensitivity factors, i.e., dRpR/dT85 min (see Table 1), 
given in Eqs. 2a, 2b, and 2c have to be optimized or tuned with the help of radar 
observations. This tuning procedure is discussed in the following section. 

3. Tuning Procedure for the Rain Retrieval Model 

For the purpose of optimization or tuning, we use a dependent sample of six 
rain events on land, and another sample of six rain events on ocean. For land, 
details pertaining to these events are listed in Tables 2a and 2b. For ocean, these 
details are listed in Tables 3a and 3b. 

In the initial stage of the tuning procedure, preliminary estimates of the 
sensitivity factors needed in the retrieval algorithm are deduced as follows. For a 
given rain event, we identify the positions and the categories of Cbs with the help of 
T85 min and |dT85/dr| (see for e.g. Fig. Id). Then, for each Cb in that event, we can 

compute Rpr from the PR data. A similar procedure is applied to the remaining 
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rain events, and a data set consisting of Cbs and their corresponding values of Rpr 
is developed. From this data set, we isolate one subset consisting of young Cbs. For 
that one subset, we determine a best linear fit between T85 m in and Rpr . The slope of 
this linear fit for that Cb category gives us a preliminary estimate of the sensitivity 
factor {dRpR/dT85 min }y. In an analogous fashion, we determine preliminary 
values of the sensitivity factors {dR PR /dT85 min } M and{dRp R /dT85 min } D applicable 
to mature and decaying Cbs, respectively. The preliminary estimates of the 
sensitivity factors are then refined as follows. 


Using Eqs. 1, 2a-2c, and 3, with the given sensitivity factors, we compute the 
spatial distribution of rain rates R for the 12 dependent mesoscale rain events 3 
mentioned earlier. Then, utilizing the estimated values of R in the 2°x3° region 
that overlap with the PR scan area, the following seven quantities are calculated for 
each rain event: 

• f^/ an d f ^ - quantities directly related to frequency of rain rate in three 
intervals, I) 1-10 mmhr'l, II) 10-20 mmhr'l, and HI) greater than or equal to 
20 mmhr'l. These three quantities are simply the fractional rain areas 

corresponding to these three rain intervals.; 

• (R)j/ (R) n , and (R) m - mean rate rates in the three intervals; 

• (R) * ■ area-average rain rate in the mesoscale region. 


Note that these seven quantities are not independent of one another. They are 
inter-dependent. 


3 These mesoscale rain events are represented using grid boxes of 2° lat. x 3° Ion.. This choice of grid box 
size is compatible with the 220 km swath of the PR. 
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To improve the preliminary values of the sensitivity factors, the ensemble 
means of each of the above seven quantities, denoted as (R)j, (R)jj, (R)jjj, f\ f H , 

f , and (R) a , are computed from the dependent sample of six rain events 

separately for the land and ocean MCS cases. These ensemble means are compared 
with those given by the PR. Based on this comparison, adjustments are made to the 
sensitivity factors in Eqs. 2a, 2b, and 2c until these seven ensemble-mean quantities 
given by the radiometer and the corresponding quantities given by the radar 
converge to within about 15 % both on land and ocean. 

We may note that the values of these sensitivity factors will depend on the 
criteria chosen for T85 m in an d | d T85/ dr | , which are used to separate the three Cb 

categories (see Fig. 3). The above procedure gives us an optimized rain retrieval 
method that can be applied to a mesoscale region. In Tables 2a, 2b, 3a, and 3b, we 
show the results obtained from the six dependent MCS cases separately over land 
and ocean. However, since this sample of rain events is used in the tuning 
procedure, we have to demonstrate the applicability of this method with an 
independent set of rain events. 

4. Application to Independent Cases 

In order to show the general applicability of our algorithm, we have retrieved 
rain rates, R, for 10 independent cases over land, and for the same number of cases 
over ocean. Also, over land and ocean, the TRMM Standard Rain Rates, Rpr and 
Rtml have been assembled for these cases. Here, the rain rates Rpr and Rtmi are the 
Version 4 Official Standard Products of the TRMM Project. The relevant details 
pertaining to these events are listed in Tables 2a and 2b for land, and in Tables 3a 
and 3b for ocean. 
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The rain retrievals, Rtmi and R/ can be scrutinized with the mesoscale area- 
average rain rates, (R^,^ and (R) A - We find (R) A and (R tm ) a given in 

Tables 2b and 3b correlate well with those given by • However, this result 

has limited significance. The correlation of mesoscale area-average rain rate does 
not guarantee quantitative agreement of the stratiform and convective rain 
distributions, and the associated latent heat release into the atmosphere (see for e.g., 
Simpson et al., 1996). 

The merits of our retrieval algorithm can be assessed by comparing the 
frequency distribution of rain rate. This comparison is done in a gross fashion in 
this study as discussed in Section 3 by considering frequency of rain rate in three 
intervals: 1-10 mmhr‘1 (f ), 10-20 mmhr‘1 (f^) and greater than or equal to 
20 mmhr' 1 (f ). From Tables 2a and 3a, the frequency of rain rate in these three 
intervals deduced from R are shown to agree to within 15 % with those given by the 
PR both on land and ocean. On land, the TMI Standard Rain Rates, RtmJ/ are not 
tuned to a radar. As a consequence, the frequency distribution given by Rtmi on 
land shows a poor agreement with PR. On ocean, the TMI Standard Rain Rates 
have been tuned to a radar. Thus, the lack of agreement in the frequency 
distribution of Rtmi over ocean reveals a weakness of the TMI Standard Rain Rates. 
This weakness has been pointed out in an independent assessment by Olson et al. 
(1999) (see Introduction). 

In Tables 2b and 3b, we note (R)jjj agrees to within 15 % of {R P r)jjj • Such a 

comparison of mean rain rates in the intervals I and II, not shown in the tables, also 
reveal an agreement to within 15 %. Similar statistics based on Rtmi show a poor 
agreement with those of the PR. From the independent sample of rain events, we 
contend that the algorithm developed here is capable of producing over land and 
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ocean mean rain rates (R)j, (R)jp ar id an d their associated frequency 

distribution, that are in better agreement with that given by the radar. 

In order to illustrate the properties of rain retrievals more vividly, we are 
presenting in Figures 4a-4c and 5a-5c maps of rain rate Rpr, Rjmi, and R for the 
following independent rain events on land: 1) an extensive MCS over Northern 
Argentina and Paraguay on 28 December 1998 (Day 362); and 2) an intense squall line 
with densely-packed convective cloud clusters over the Southern United States on 
11 June 1998 (Day 162). As remarked earlier, since Rtmi is not tuned to radar on 
land, it does not reveal the convective features as seen in the maps of Rpr. On the 
other hand, the rain maps of R derived from our scheme show the convective 
features in better agreement with Rpr. In order to be comprehensive, we have 
deduced a detailed frequency distribution, or histogram, of rain rate from the maps 
of Rpr, Rtmi, and R. We have also deduced for each rain event f*, f^, and f^. 
Histograms, and numerical values of f^, f^, and , for the two rain events on 
land are shown in Figures 4d-4f and 5d-5f. 

We may point out that the rain map shown in Figure 5c, derived from our 
scheme, reveals reasonably well the intense rain from densely-packed convective 
cloud clusters observed along a squall line over the Southern United States on 
11 June 1998 (Day 162). Some closely-spaced Cbs in Figure 5c form convective areas 
that have dimensions of about 40 km. The geographic distribution of intense rain 
given by our algorithm agrees well with that given by Rpr. 

To illustrate the ability of our rain retrieval scheme over ocean, in Figures 6a- 
6c and 7a-7c we are presenting maps of rain rate Rpr, Rtmi, and R for two cases: 1) 


18 


LATITUDE LATITUDE LATITUDE 





LONGITUDE 



LONGITUDE 




D 
O 

Lil 

£ 20 * 


f =70 % 
f=18 % 
f =8 % 

<R> a = 10.1 mmhr 1 


lot 




20 


AO 


100 


R mmhr 1 


Figure 4: Maps of surface rain rate, and its frequency distribution, corresponding to a rain event that 
occurred over Northern Argentina and Paraguay on 28 December 1998 (Day 362). 

a) map of PR rain rate, RpR (mmhr-1). 

b) map of TRMM Standard Rain Rate, RyMI (mmhr-1), deduced from the microwave radiometer. 

c) map of rain rate, R (mmhr-1), retrieved from our microwave radiometer algorithm. 

d) frequeny distribution of PR rain rate, RpR (mmhr-1). 

e) frequeny distribution of TRMM Standard Rain Rate, RtMI (mmhr-1), deduced from the microwave 
radiometer. 

f) frequeny distribution of rain rate, R (mmhr-1), retrieved from our microwave radiometer algorithm. 
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Figure 5: Maps of surface rain rate, and its frequency distribution, corresponding to a rain event that 
occurred over the Southern United States on June 11, 1998 (Day 162). 

a) map of PR rain rate, Rj>r (mmhr-1). 

b) map of TRMM Standard Rain Rate, RjMI (mmhr-1), deduced from the microwave radiometer. 

c) map of rain rate, R (mmhr-1), retrieved from our microwave radiometer algorithm. 

d) frequeny distribution of PR rain rate, RpR (mmhr-1). 

e) frequeny distribution of TRMM Standard Rain Rate, Rx\n(mmhr-1), deduced from the microwave 
radiometer. 

f) frequeny distribution of rain rate, R (mmhr-1), retrieved from our microwave radiometer algorithm. 
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Figure 6: Maps of surface rain rate, and its frequency distribution, corresponding to a rain event that 
occurred over the Western Equatorial Pacific on January 24, 1999 (Day 24). 

a) map of PR rain rate, RpR (mmhr-1). 

b) map of TRMM Standard Rain Rate, RtmI (mmhr-1), deduced from the microwave radiometer. 

c) map of rain rate, R (mmhr-1), retrieved from our microwave radiometer algorithm. 

d) frequeny distribution of PR rain rate, RpR (mmhr-1). 

e) frequeny distribution of TRMM Standard Rain Rate, RTMl( mm h r_ l)> deduced from the microwave 
radiometer. 

f) frequeny distribution of rain rate, R (mmhr-1), retrieved from our microwave radiometer algorithm. 







Figure 7: Maps of surface rain rate, and its frequency distribution, corresponding to Hurricane Floyd 
as it moved across the Western Atlantic on September 13, 1999 (Day 256). 

a) map of PR rain rate, RpR (mmhr-1). 

b) map of TRMM Standard Rain Rate, RpMI (mmhr-1), deduced from the microwave radiometer. 

c) map of rain rate, R (mmhr-1), retrieved from our microwave radiometer algorithm. 

d) frequeny distribution of PR rain rate, RpR (mmhr-1). 

e) frequeny distribution of TRMM Standard Rain Rate, RiMl(mmhr-l), deduced from the microwave 
radiometer. 

f) frequeny distribution of rain rate, R (mmhr-1), retrieved from our microwave radiometer algorithm. 








the 24 January 1999 (Day 24) rain event over the Equatorial Western Pacific; and 2) 
Hurricane Floyd as it moved across the Western Atlantic Ocean on 13 September 
1999 (Day 256). The Hurricane Floyd case is not one of the cases presented in 
Tables 3a and 3b. The histograms and numerical values f*, and of rain rate 
given by Rpr, Rtml and R are shown in Figures 6d-6f and 7d-7f, respectively. From 
an examination of these figures we find the Standard Rain Rate R TMI , when 
compared to Rpr, is not as good as R on the ocean. 

The rain retrieval algorithm developed in this study is tuned to MCS cases 
over land and ocean. Such tuning does not guarantee its applicability to a hurricane 
event. The rain map of Hurricane Floyd shown in Fig. 7c, derived from our 
algorithm, captures in essence all the gross features found in the rain map of Rpr 
(Fig. 7a). For example, contrasts between large rain rates in the eyewall and spiral 
bands and light rain rates in the weak stratiform areas are captured better in our 
algorithm. Also, the histograms of Rpr (Fig. 7d) and R (Fig. 7f) have a log-normal 
character, while that given by Rtmi (Fig- 7e) does not. 

In order to obtain a satisfactory convective and stratiform rain distribution 
from the measurements made by microwave radiometers, we find empirical tuning 
with collocated radar rain observations is needed. For this reason, the empirical 
rain retrieval scheme of Hong et al. (1999) has been tuned to a radar. This scheme is 
adopted to retrieve TMI Standard Rain Rates on ocean. However, Olson et al. (1999) 
point out the weakness of this algorithm in capturing quantitatively convective and 
stratiform rain over ocean. The reason for this weakness is that the scheme of Hong 
et al. (1999) does not have a way to limit convective rain area in an MCS to a 
thunderstorm scale. In the empirical rain retrieval algorithm developed in this 
study, convective rain in Cbs is limited to a thunderstorm scale. Hong et al. have 
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not demonstrated the utility of their algorithm by tuning it with radar on land. Our 
rain retrieval scheme, developed with the help of the PR, can satisfactorily produce 
convective and stratiform rain rate distributions over both land and ocean. 

The partitioning of convectively active thunderstorms, i.e., young and 
mature Cbs, from those that are decaying is critical for the determination of latent 
heat input into the atmosphere. In Figure 8, for all the MCS rain events over land 
and ocean, we show separately the frequency distributions of these active and 
decaying Cbs as a function of T85 m in- From this figure, we see that the frequency 
distributions of the convectively active Cbs, both on land and ocean, tend to peak at 
a much colder T85 m i n than the decaying Cbs. This vividly illustrates the stronger 
scattering effect of the dense ice particles in the active Cbs versus the less-dense ice 
a gg re g ates in the decaying Cbs. This simple analysis presented in Fig. 8 lends 
support to the criterion jdT85/dr| used to classify the convectively active and 
decaying Cbs (see Fig. 3). 

The present algorithm has the ability to retrieve rain in relatively warm 
clouds cells of young Cbs that have sufficient vertical growth to produce some ice. 
Purely warm rain where there is no ice formation will be missed in our algorithm. 
By comparing the PR rain with the rain retrieved from our algorithm, such missing 
warm rain can be identified. From such a comparison from several MCS events, we 
find missing warm rain area is less than about 5 % of the total rain area. 

Based on the results and the comparison with the radar observations 
presented above, we note that the radiometer derived rain information can capture 
crudely the convective and stratiform rain character given by the PR. However, it is 
not capable of resolving fine details. The main advantage of our TMI rain retrieval 
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Figure 8: Frequency distribution of the convectively active young and mature thunderstorms, as well as the decaying thunderstorms, 
over land and ocean. 



algorithm is that it can extend in a crude way the PR quality rain information to the 
760 km swath width of the radiometer. As a result, one can obtain better estimates 
of latent heat release into the atmosphere over a wider area, which in turn can be 
used to improve general circulation models. 

5. Conclusions 

Multi-spectral, dual-polarization measurements made by satellite-borne 
microwave radiometers do not contain enough independent pieces of information 
to derive the amount and the vertical distribution of liquid, mixed phase, and 
frozen hydrometeors present in the atmosphere. This leads to uncertainties in 
quantitative estimation of rain rate pertaining to convective and stratiform rain 
based purely on radiative transfer theoretical considerations. Additional 
information, that is independent of the spectral measurements, is necessary to 
derive rain quantitatively. In Prabhakara et al. (1998, 1999), the fractional rain area 
in a mesoscale region derived from the microwave radiometer is used as additional 
information to derive mesoscal e-average rain rate. The basis for this broad 
relationship is contained in the dynamics of Mesoscale Convective Systems which 
produce rain in that area. This method, based on fractional rain area, does not 
identify stratiform from convective rain areas. 

In the last two decades, remote sensing of rain with the help of satellite-borne 
microwave radiometer observations has made significant progress. At present, the 
spatial resolution of the Tropical Rainfall Measuring Mission (TRMM) Microwave 
Imager (TMI) radiometer (~ 5.5 km at 85 GHz) has improved to the point that it can 
readily reveal some detailed information on a scale of thunderstorms, i.e., about 
20 km. This capability is needed to retrieve convective rain. In the present method. 
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we utilize TMI radiometer data, and the conceptual framework that thunderstorm 
scale dynamics (~ 20 km) establishes a relationship between the spatial distribution 
of T85 and the existence of Cbs of different types: young, mature , and decaying . A 
significant result of this study is that we can discriminate convectively active young 
and mature Cbs from stratiform or decaying Cbs. Also, we take into consideration 
radiative transfer properties (see for e.g. Schols et al., 1999) associated with the 
hydrometeors of these three different categories of Cbs, and infer relationships 
between the rain rate in these Cbs and T85 m i n . Initially, the rain rate observations 
made by the TRMM PR are used to tune these relationships. This entire procedure 
constitutes the basis for separation of convective rain from stratiform rain. 

Results produced from independent rain events over land and ocean show a 
good consistency between our method and the PR. As a consequence of this 
capability, we can extend the rain information given by the 220 km swath of the PR 
to the full 760 km swath of the TMI. The method developed here is 
computationally simple and fast, and can be used as an effective tool to generate the 
climatology of precipitation over tropical land and ocean. 

We may bring to the attention of the reader some assumptions and 
weaknesses in the present rain retrieval method. We are assuming that all Cbs, 
regardless of their type, have a circular shape. Hie strong updrafts in the Cbs may be 
tilted, which leads to a given T85 minimum that does not collocate with the actual 
center of the Cb. This mislocation is ignored. The actively developing and decaying 
portions of a Cb that may be close together (Houze, 1997) are not resolved in the 
present method. At the wavelength used by the PR (~ 2 cm), attenuation at high 
rain rates can be significant. When the rain rate is high, this can introduce errors in 
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the PR rain rate estimates. Since our algorithm is timed to the PR data, errors in the 
PR rain estimates will be reflected in our retrievals. 

The empirical rain retrieval method developed here, with the aid of the PR, 
can serve as a useful tool for retrieving rain over the 760 km swath of the TMI. 
However, this radiometer derived rain information can capture only the gross 
character of the radar observations, and not the fine details. Improvement in the 
spatial resolution of the microwave radiometer, and better methods to identify 
different kinds of hydrometeors, will be helpful in advancing the remote sensing 
capability. 
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Figure Captions 


Figure 1: Data pertaining to a Mesoscale Convective System over the Southeast 
United States on June 5, 1998 (Day 156). 

a) Map of 85 GHz brightness temperature, T85 (K), deduced from TMI. 

b) Map of TRMM Precipitation Radar surface rain rate Rpr (mrnhr" 1 ). 

c) Frequency distribution of surface rain rate for this event deduced from PR data. 

d) Map of minima of T85 for this rain event indicating location of different types of 
thunderstorms - young (Y), mature (M) and decaying (D). The numbers shown on 

the left and right of the symbols Y, M and D indicate the magnitude of 
| dT85/dr | (Kknr 1 ) and T85 m in (K). 

Figure 2: Illustration of the method to determine the location of a T85 minimum, 
T85 m in, and the mean spatial gradient, |dT85/dr|, around this minimum. The solid 

curves represent the conical scan of the TMI radiometer, and the circles represent 
the footprints. Along the scan, the separation between footprints is about 4.6 km, 
while the separation between scans is about 13.9 km. In the illustration, the 
coefficients A and B correspond to appropriate weighting factors. 


Figure 3: Scatterplot of |dT85/dr| versus T85 m i n corresponding to decaying (D), 
young (Y), and mature (M) Cbs. The horizontal line ab, where |dT85/dr| is equal 

to 1 Kkm'l, separates the D Cbs from those of the Y and M categories. The vertical 
line ef , where T85 min is equal to 210 K, separates the Y from the M Cbs. 



Figure Captions Continued 


Figure 4: Maps of surface rain rate, and its frequency distribution, corresponding to a 
rain event that occurred over Northern Argentina and Paraguay on 28 December 
1998 (Day 362). 

a) map of PR rain rate, Rpr (mmhr'l). 

b) map of TRMM Standard Rain Rate, Rjmi (mmhr‘1), deduced from the microwave 
radiometer. 

c) map of rain rate, R (mmhr 1 ), retrieved from our microwave radiometer 
algorithm. 

d) frequency distribution of PR rain rate, Rpr (mmhr' 1 ). 

e) frequency distribution of TRMM Standard Rain Rate, Rtmi (mmhr'l), deduced 
from the microwave radiometer. 

f) frequency distribution of rain rate, R (mmhr**), retrieved from our microwave 
radiometer algorithm. 


Figure Captions Continued 


Figure 5: Maps of surface rain rate, and its frequency distribution, corresponding to a 
rain event that occurred over the Southern United States on June 11, 1998 (Day 162). 

a) map of PR rain rate, Rpr (mmhr 1 ). 

b) map of TRMM Standard Rain Rate, Rtmi (mmhr' 1 ), deduced from the microwave 
radiometer. 

c) map of rain rate, R (mmhr -1 ), retrieved from our microwave radiometer 
algorithm. 

d) frequency distribution of PR rain rate, Rpr (mmhr- 1 ). 

e) frequency distribution of TRMM Standard Rain Rate, Rtmi (mmhr* 1 ), deduced 
from the microwave radiometer. 

f) frequency distribution of rain rate, R (mmhr* 1 ), retrieved from our microwave 
radiometer algorithm. 



Figure Captions Continued 


Figure 6: Maps of surface rain rate, and its frequency distribution, corresponding to a 
rain event that occurred over the Western Equatorial Pacific on January 24, 1999 
(Day 24). 

a) map of PR rain rate, Rpr (mmhr 1 ). 

b) map of TRMM Standard Rain Rate, Rtmi (mmhr'l), deduced from the microwave 
radiometer. 

c) map of rain rate, R (mmhr'l), retrieved from our microwave radiometer 
algorithm. 

d) frequency distribution of PR rain rate, Rpr (mmhr'l). 

e) frequency distribution of TRMM Standard Rain Rate, Rtmi (mmhr' 1 ), deduced 
from the microwave radiometer. 

f) irequency distribution of rain rate, R (mmhr'l), retrieved from our microwave 
radiometer algorithm. 


Figure Captions Continued 


Figure 7: Maps of surface rain rate, and its frequency distribution, corresponding to 
Hurricane Floyd as it moved across the Western Atlantic Ocean on September 13, 
1999 (Day 256). 

a) map of PR rain rate, Rpr (mrnhr'l). 

b) map of TRMM Standard Rain Rate, Rtmi (mmhr'l), deduced from the microwave 
radiometer. 

c) map of rain rate, R (mmhr 1 ), retrieved from our microwave radiometer 
algorithm. 

d) frequency distribution of PR rain rate, Rpr (mmhr'l). 

e) frequency distribution of TRMM Standard Rain Rate, R^ (mmhr' 1 ), deduced 
from the microwave radiometer. 

f) frequency distribution of rain rate, R (mmhr'l), retrieved from our microwave 
radiometer algorithm. 

Figure 8: Frequency distribution of the convectively active young and mature 
thunderstorms, as well as the decaying thunderstorms, over land and ocean. 



Table 1: Sensitivity of T85 m i n to the Cb mean surface rain rate Rpr 
tropical land and ocean regions. 

Sensitivity Factors 

CbType 85 GHz Minimum dRpR/dT85 min 


Young 255 K > TSSmin > 210 K ' 0.25 mmhr' 1 K' 1 

Mature T85nun < 210 K 0.35 mmhr'lK‘1 

Decaying 255 K > T85 min 0.12 miuhrlK* 1 


applicable to 



Table 2a: Land dependent and independent MCS events and the associated fractional rain area 


statistics. The fractional rain area is calculated for each of the following rain rate intervals: I) 1-10 


mmhr"!. 

II) 10-20 mmhr'l. 

and III) greater than or 

equal to 

20 mmhr'l. 
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0.21 
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0.03 

0.02 

0.0 

169/98 

26.9 

115.7 

0.57 

0.51 

0.35 

0.04 

0.03 

0.25 

0.02 

0.01 

0.0 

171/98 

25.8 

115.0 

0.22 

0.19 

0.17 

0.01 

0.02 

0.02 

0.02 

0.01 

0.0 

180/98 

32.8 

160.0 

0.27 

0.28 

0.11 

0.04 

0.06 

0.21 

0.02 

0.04 

0.0 

201/98 

28.9 

116.0 

0.29 

0.25 

0.07 

0.02 

0.02 

0.18 

0.03 

0.01 

0.0 

MEAN 



0 . 32 

0.27 

0.14 

0.03 

0.03 

0 . 17 

0,02 

0 . 02 

0 . 00 






Independent Cases 






Day/Yr 

Lat 

Lon 

f 1 

f 1 

r PR 

f 

J TMI 

f H 

f 11 

*PR 

f H 

a tmi 

£ III 

r m 
r PR 

f m 

a tmi 

146/98 

35.0 

- 89.5 

0.48 

0.38 

0.17 

0.06 

0.11 

0.37 

* 

0.04 

0.08 

0.0 

156/98 

35.0 

- 89.1 

0.62 

0.47 

0.12 

0.10 

0.13 

0.56 

0.06 

0.09 

0.02 

160/98 

4.4 

- 66.2 

0.45 

0.47 

0.39 

0.02 

0.02 

0.05 

0.0 

0.0 

0.0 

162/98 

33.0 

- 99.5 

0.29 

0.19 

0.15 

0.08 

0.02 

0.19 

0.06 

0.09 

0.06 

162/98 

31.7 

115.6 

0.36 

0.31 

0.25 

0.02 

0.01 

0.04 

0.0 

0.0 

0.0 

186/98 

24.0 

81.5 

0.44 

0.49 

0.18 

0.06 

0.06 

0.31 

0.03 

0.03 

0.0 

271/98 

31.5 

- 87.0 

0.76 

0.72 

0.48 

0.06 

0.10 

0.35 

0.02 

0.03 

0.00 

362/98 

- 26.7 

- 58.0 

0.70 

0.52 

0.06 

0.18 

0.16 

0.86 

0.08 

0.12 

0.03 

19/99 

- 17.9 

29.8 

0.32 

0.40 

0.19 

0.04 

0.03 

0.12 

0.03 

0.02 

0.0 

44/99 

- 17.2 

127.5 

0.43 

0.28 

0.15 

0.05 

0.05 

0.31 

0.06 

0.04 

0.01 

MEAN 



0.48 

0.42 

0.21 

0.07 

0.07 

0.32 

0 . 04 

0 .05 

0 .01 


I 


Table 2b: Land dependent and independent MCS events and the associated mesoscale-mean rain 
statistics. The mean rain rate is calculated for each of the following rain rate intervals: 

I) 1-10 mmhr'l, II) 10-20 mmhr'l, and III) greater than or equal to 20 mmhr _ L 


Dependent Cases 


Day/Yr Lat Lon 


156/98 32.8 - 92.0 
159/98 12.5 - 1.4 
169/98 26.9 115.7 
171/98 25.8 115.0 
180/98 32.8 160.0 
201/98 28.9 116.0 

MEAN 


Win 

( RpR )lII 

(^TMl)j 

31.3 

32.0 

0.0 

37.3 

30.5 

20.0 

33.4 

30.0 

0.0 

35.7 

34.9 

0.0 

30.3 

30.5 

0.0 

32.8 

25.1 

0.0 

33.5 

30.5 

3.3 


( R pr) a ( R tmi) a 

1.7 1.9 2.0 

3.2 1.4 3.2 

2.9 2.3 5.1 

1.3 1.1 1.2 

2.1 3.1 3.0 

2.1 1.1 2.6 

2.2 1.8 2.9 


Independent Cases 


Day/Yr 

Lat 

Lon 

(»)„, 

W,n 

( R TMl) ln 

w » 

( r pr) a 

(^Ml) 

146/98 

35.0 

- 89.5 

39.4 

38.2 

0.0 

4.4 

6.3 

5 . 6 

156/98 

35.0 

- 89.1 

40.2 

44.2 

25.5 

6.8 

7.7 

8 4 

160/98 

4.4 

- 66.2 

26.0 

32.3 

0.0 

1.5 

1.7 


162/98 

33.0 

- 99.5 

61.6 

75.1 

20.8 

5.7 

7.8 

z . o 

4 . 8 

162/98 

31.7 

115.6 

24.8 

21.5 

0.0 

1.2 

,7 

n 

186/98 

24.0 

81.5 

33.3 

28.4 

0.0 

3.4 

3 .4 

Z . U 

4 8 

271/98 

31.5 

- 87.0 

29.3 

27.0 

0.0 

3.6 

5.1 

1 ■ o 

7 1 

362/98 

- 26.7 

- 58.0 

42.1 

45.2 

20.9 

10.1 

9.9 

12 . 1 

19/99 

- 17.9 

29.8 

36.8 

28.8 

20.0 

2.4 

2.0 

2 6 

44/99 

- 17.2 

127.5 

42.3 

32.8 

20.5 

5.0 

3.0 

5.2 

MEAN 



37 . 6 

37.4 

10.8 

4.4 

4.8 

5 . 5 


Table 3a: Ocean dependent and independent MCS events and the associated fractional rain area 
statistics. The fractional rain area is calculated for each of the following rain rate intervals: 

I) 1-10 mmhr'l, II) 10-20 mmhr 1 , and III) greater than or equal to 20 mmhr 1 . 


Dependent Cases 


Day/Yr 

Lat 

Lon 

f 1 

fL, 

f 1 


f 11 

f 11 

f m 

rin 

f m 





1 PR 

A TMI 


r PR 

i TMI 

I 

r PR 

1 TMI 

159/98 

- 9.5 

85.0 

0.33 

0.44 

0.46 

0.02 

0.02 

0.17 

0.01 

0.00 

0.01 

162/98 

13.6 

- 105.2 

0.41 

0.39 

0.51 

0.02 

0.02 

0.06 

0.01 

0.01 

0 . 03 

184/98 

3.9 

- 91.8 

0.47 

0.45 

0.56 

0.04 

0.06 

0.15 

0.03 

0.03 

0 . 05 

314/98 

5.5 

147.2 

0.63 

0.61 

0.62 

0.06 

0.08 

0.35 

0.04 

0.03 

0.06 

7/99 

4.5 

156.0 

0.45 

0.46 

0.55 

0.04 

0.03 

0.14 

0.01 

0.01 

0.03 

70/99 

- 25.2 

54.0 

0.59 

0.58 

0.54 

0.02 

0.02 

0.20 

0.0 

0.00 

0.01 

MEAN 



0.48 

0 . 49 

0.54 

0.03 

0.04 

0.18 

0 . 02 

0 .01 

0 .03 


Independent Cases 


Day/Yr 

Lat 

Lon 

f 1 

f* 

f 1 





a pr 

a tmi 

154/98 

6.3 

- 13.0 

0.59 

0.47 

0.64 

158/98 

21.0 

118.2 

0.38 

0.28 

0.48 

244/98 

26.1 

- 90.6 

0.45 

0.40 

0.59 

269/98 

8.5 

- 129.8 

0.41 

0.36 

0.46 

335/98 

1.0 

147.0 

0.31 

0.24 

0.34 

24/99 

6.0 

148.5 

0.64 

0.65 

0.43 

63/99 

- 0.4 

- 36.0 

0.55 

0.36 

0.45 

65/99 

- 14.4 

102.4 

0.71 

0.64 

0.62 

68/99 

- 27.0 

- 29.6 

0.38 

0.27 

0.47 

73/99 

- 25.9 

168.8 

0.59 

0.47 

0.54 

MEAN 



0 .50 

0.41 

0 .50 


f H 

f 11 


fill 

fin 

f m 


X PR 

1 TMI 

1 

r PR 

a tmi 

0.02 

0.04 

0.25 

0.01 

0.02 

0.01 

0.01 

0.01 

0.04 

0.01 

0.02 

0.02 

0.04 

0.06 

0.14 

0.02 

0.02 

0.04 

0.03 

0.03 

0.13 

0.0 

0.0 

0.01 

0.01 

0.01 

0.03 

0.01 

0.01 

0.01 

0.08 

0.06 

0.38 

0.05 

0.05 

0.08 

0.05 

0.06 

0.14 

0.02 

0.02 

0.03 

0.06 

0.05 

0.12 

0.04 

0.06 

0.12 

0.06 

0.05 

0.13 

0.03 

0.04 

0.06 

0.01 

0.01 

0.06 

0.01 

0.01 

0.02 

0.04 

0 . 04 

0.14 

0 . 02 

0 . 02 

0 . 04 




Table 3b: 

Ocean dependent and independent MCS events and the associated 

mesoscale- 

■mean 

rain statistics. The mean rain rate is calculated for each of the following rain rate intervals: 

I) 1-10 mmhr 1 , 

II) 10-20 mmhr'i, 

and III) greater than or equal to 20 mmhr“l 






Dependent 

Cases 




Day/Yr 

Lat 

Lon 

< r >iii 

( RpR ) III 

( r tmi) iii 

«» 

( RpR ) A 

( r tmi) a 

159/98 

- 9.5 

85.0 

21.8 

20.0 

24.5 

1.5 

1.6 

4.2 

162/98 

13.6 

- 105.2 

30.2 

28.5 

26.9 

1.9 

1.5 

3.3 

184/98 

3.9 

- 91.8 

34.4 

29.5 

31.1 

3.3 

3.2 

6.0 

314/98 

5.5 

147.2 

30.5 

28.6 

25.1 

4.2 

4.4 

8.3 

7/99 

4.5 

156.0 

28.2 

32.4 

22.7 

2.3 

2.4 

4.5 

70/99 

- 25.2 

54.0 

20.3 

21.2 

24.6 

2.3 

1.9 

4.9 

MEAN 


4-‘ 

27 . 6 

26.7 

25.8 

2 . 6 

2 . 5 

5 . 2 


i 





Independent 

Cases 




Day/Yr 

Lat 

Lon 

Win 

( r pr)hi 



( r pr) a 

( r tmi) a 

154/98 

6.3 

- 13.0 

29.3 

29.4 

25.3 

2.0 

2.5 

5.7 

158/98 

21.0 

118.2 

32.5 

37.9 

25.8 

1.5 

1.6 

2.2 

244/98 

26.1 

- 90.6 

30.8 

30.4 

27.3 

2.6 

2.8 

5.1 

269/98 

8.5 

- 129.8 

22.3 

22.9 

21.8 

1.7 

1.6 

3.6 

335/98 

1.0 

147.0 

28.8 

37.4 

23.9 

1.2 

1.0 

1.7 

24/99 

6.0 

148.5 

33.0 

33.0 

25.4 

5.3 

5.1 

9.6 

63/99 

- 0.4 

- 36.0 

30.3 

25.2 

26.2 

3.4 

2.7 

4.7 

65/99 

- 14.4 

102.4 

40.6 

44.2 

40.1 

4.9 

4.9 

8.7 

68/99 

- 27.0 

- 29.6 

29.4 

35.0 

24.9 

3.2 

3.0 

5.0 

73/99 

- 25.9 

168.8 

37.5 

28.3 

28.8 

1.9 

1.6 

2.6 

MEAN 



31.4 

32.4 

27.0 

2 . 8 

2.7 

4 . 9 



LONGITUDE 





LONGITUDE 




o 

o 



o 

00 


o 

to 



o 

CM 


o 


I c. 


mmhr 1 








When T 0 <T 1 ,T 2 , T 3 and T 4 
Then T m j n =TQ 

A(T-| +T 2 ) + B(T 3 + T 4 ) - 4Tq 

5- a 


z 



3 


min 















LONGITUDE 




25.60 



aanuivn 


LONGITUDE 



25.60 



aamiivn 


LONGITUDE 








R mmhr 1 














mmhr 










aanmvn 






mmhr 1 







mmhr 1 


















mmhr 1 










